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Abstract: In this note, we establish an asymptotic expansion for the cen- 
tering parameter appearing in the central limit theorems for linear spectral 
statistic of large-dimensional sample covariance matrices when the pop- 
^ ' ulation has a spiked covariance structure. As an application, we provide 

\| ' an asymptotic p ower function for the corrected likelihood ratio test in 

>^3 ' iBai et ahl l|200y) used for testing the presence of spike eigenvalues in the 

population covariance matrix. This result generalizes a formula provided in 
VP , lOnatski et al. 1 1I2OIII ) where only one simple spike exists. 

\l . AMS 2000 subject classifications: Primary 60F05; secondary 62H15. 

(a_>^ ' Keywrords and phrases: Large-dimensional sample covariance matrices, 

CO ' Spiked population model, Central limit theorem. Centering parameter, fac- 

tor models. 



1. Introduction 

Let (Sp) be a sequence oi px p non-random and nonnegative definite Hermi- 
tian matrices and let (wij), i, j > 1 be a doubly infinite array of i.i.d. complex- 
valued random variables satisfying 

E(«;u) = 0, E(|«;ii|2) = l, E{\wnf)<^. 

Write Zn = (wij)i<i<p^i<j<„, the upper-left p x n block, where p = p{n) is 
related to n such that when n — >■ (X), p/n — > y > 0. Then the matrix Sn = 

—Yip Z„Z*^'Sp can be considered as the sample covariance matrix of an i.i.d. 

1 /? 
sample (xi,...,x„) of p-dimensional observation vectors Xj = Ep Uj where 
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Uj = {'Wij)i<i<p denotes the j-th column of Z„. Note that for any semi-positive 
definite Hermitian matrix A, A^^^ denotes a Hermitian square root and we cah 
the spectral distribution (SD) the distribution generated by its eigenvalues. 

Assume that the SD iJ„ of Sp converges weakly to a nonrandom probabil- 
ity distribution H on [0, oo). It is then well-known that the SD i^'^" of 5'„, 
generated b y its eigenvalues An.i > • • ■ > A^^p, converges t o a nonrandom lim- 



iting SD G (jMarcenko and Pastuii Il967t ISilversteinl 119951) . The so-called null 
case corresponds to the situation Ep = /p, so iJ„ = Si and the limiting SD is 
the seminal Marcenko-Pastur law Gy with index y and support [ay, by] where 
1 — ^/?y)^, 6„ = (1 -I- \f yf' , and an additional mass at the origin if y > 1 



(JMarcenko and Pastuj (|l967l )) 



In t his paper we consider the spiked population model introduced in lJohnstone 



(|2001l ) where the eigenvalues of Sp are 



ai, 



,"1, 



(1.1) 



Here M and the multiplicity numbers (n^) are fixed and satisfy ni + ■ • • + Wfe = 
M . In other words, all the population eigenvalues are unit except some fixed 
number of them (the spikes) . The model can be viewed as a finite-rank pertur- 
bation of the null case. Obviously, the limiting SD G of S'„ is not affected by 
this perturbation. However, the asymptotic behaviour of the extreme eigenval- 
ues of Sn is significantly different from the null case. The analysis of this new 
beh aviour of ex t reme e igenvalues has been an activ e area i n the l ast years, see 
e.g . iBaik et al.l (l2005l). Baik and Silversteinl (l2006l). IPaull dJOOTh. iBai and Ya( ' 



(l2008l).lBenavch-Georges et al. l(|201ll ). lBenavch-Georges and Nadakuditil(|201l[ 
and iBai and Yad ( 20121 ). In particular, the base component of the population 
SD Hn in the last two references has been extended to a form more general than 
the simple Dirac mass 5i of the null case. 

For statistical applications, besid es the principal c omponents analysis which 
is indeed the origin of spiked models (jjohnstond (J200l[ )). large-dimensional strict 
factor models are equivalent to a spiked population model and can be analyzed 
using the above-men tioned results. Related recent contri b utions in the a r ea in- 
clude , among others. iKritchman and Nadle^ (l2008l l2009l) lOnatskil (|2009|. I2OIOI 



20121 ). iNadakuditi and Silversteinl (J2010I ) and IPassemier and Yao I (|2012l ) and 



almost all of them concern the problem of estimation and testing the number 
of factors (or spikes) . 

In this note, we analyze the effects caused by the spike eigenvalues on the 
fluctuations of linear spectral statistics of the form 



T«(/) = E/(A„,0 = F^"(/), 



(1.2) 



i=l 



where / is a given function. Similarly to the convergence of the SD's, the pres- 
ence of the spikes does not prevent a central limit theorem for T„(/); however 
as we will see, the centering term in the CLT will be modified according to the 
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values of the spikes. As this term has no expheit form, our main resuh is an 
asymptotic expansion presented in Section [5] Section |2] explains in detail the 
contour that appears in the main results. To illustrate the importance of such 
expansions, we present in Section |4] an application for the determination of the 
power function for testing the absence versus presence of spikes. The Appendix 
contains some technical derivations. 



2. Centering parameters in the CLT of the LSS from a spiked model 

Fluctuations of linear spectral s tatistics of form (|1.2I) arc indeed covered by 



a central limit the ory developed i n jBai and Silverstein (2004). The theory was 
later improved bv |Pan and Zhoul ( 20081 ) where the restriction £'(|u'ii|^) ~ 3 



matching the real Gaussian case was removed. 

Let /i , . . . , /i be L functions analytic on an open domain of the complex plan 
including the support of the limiting SD. These central limit theorems state that 
the random vector 

(X„(/i),...,X„(/i)), 

where 

Xn{f)=p[F'"if)-Fy-'"Hf)] =pJfix)d{F^--Fy-'"-){x) , 
converges weakly to a Gaussian vector 

{Xfl : ■ ■ ■ > XfL ) 

with known mean function i?[^/] and covariance function C'ov{Xf,Xg) that 
can be calculated from contour integrals involving parameters mi^) ^'Hd H, 
where rn{z) is the companion Stieltjes transform corresponding to the limiting 
SD of 5„ = -Z*I]pZ„. If the population has a spiked covariance structure, we 
know that the limit H and rn{z) remain the same as the non-spiked case, so the 
limiting parameters -E[X/] and Cov{Xf,Xg) are also unchanged. 

It is remarked that the centering parameter ppy^'^^^f) depends on a par- 
ticular distribution F^"'^"^ which is a finite- horizon proxy for the limiting SD 
of Sji- The difficulty is that F^"'^" has no explicit form; it is indeed implic- 
itly defined through rriniz) (the finite counterpart oi m(z)), which solves the 
equation: 

z = -— + ynf -^ — dHnit) . (2.3) 

This distribution depends on the SD iJ„ which in turn depends on the spike 
eigenvalues. 

More precisely, the SD iJ„ of Sp is 



H. . ^-^S 



k 

i + iyn,5a, . (2.4) 



P P ,:=l 
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The term 



k 

P 



1 

- ^ n.iSai (t) 



i=l 



vanishes when p tends to infinity, so it has no influenee when considering hmiting 
spectral distributions. However for the CLT, the term pFi'"'-f^"(/) has a p in 
front, and - X]i=i iT-i^at times p is of order 0(1), thus cannot be neglected. 



k 

i— 

Our main result is an asymptotic expansion for this centering parameter. 



Theorem 1. (large spikes) Suppose the population has a spiked population 
structure as stated in ()1.1|) . where the spike eigenvalues a[s satisfy \ai — l\ > y/y. 
For any f analytic on an open domain including the support of the Marcenko- 



Pastur distribution Gy, we have 



py^-^-if) 



k 



1 / /(_! + -y^)[^ - y ^^^d^)dm (2.5) 

^-^wJci ™ l + m'^y^m ^ {I + a,mY ' ^ ' 

k 

1 f /f 1 , Vn xsr^ {l-ai)ni 1 y„m 



2mp Jc^ m I + m' '^ {I + a.,m){l + m)^ m {I + mY ' ^ ' 

-(1 - -)G,„(/) + - E-'/K + ^) + O(^) . (2.7) 



i=l 



Pere Gy^ (/) is the integral of f with respect to the Marcenko-Pastur distribution 
with index yn = p/n and Ci is a contour, when restricted to the real axes, 
enclosing the interval [j^,j^] (or [^^1^, -i_) y(_i_, _^] ) ^hen < 

y <1 (ory> 1). 

Theorem 2. (small spikes) When the spike eigenvalues a[s satisfy ja^ — 1| < y^, 
we have 

py-'^'-if) 



ff± + .Jt^)(J)(L_ _ y '^'"»"' )dm 



1 I r' I 1 , Vn -.sr^ {l~a.i)ni 1 y„m 



/(_± + -^^)y , ^^~:;^"' (- - -^^^)dm 

m 1 + TO ■^^ (1 + a,TO,)( 1 + TO) m (1 + to)^ 



2'Kip Jq_^ m 1 + TO ^^ (i + aim)(l + m) m (1 + to) 

+a-y)G.„(/)+o(^), 

here the Gy^ (/) and the contour Ci are the same as claimed in Theorem\^ above. 

Proof, (proof of Theorem[T]) Recall that Gy^ (/) ~ J f{x)dGy^ (x) when no spike 
exists, where Gy^ is the M — P distribution with index y„- And by the Cauchy 
integral formula, it can be expressed as ^^ § f{z)m{z)dz, where the integral 
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contour is chosen to enclose the support of Gj,„ and it's limit Gy. Besides, m(the 
companion Stieltjes transform of m) satisfies the equation; 

1 , Vn 

z = — - 



m 1 + 777, 

taking derivatives on both sides with respect to z, we get: 

dz = ( — - - — — ^^^—-:r)dm . 
^m? (1 + m)2' — 

Combine these up and use the relationship between m and m. wc have: 
GyM) = -"^^J fiz)m{z)dz 

p2TTiJc, m l + m'—^ '^rn? {I + m)^ ' — ' 

where the second equality is due to the change of variable from z to m. As a 
result, the integral region, which encloses the support of Gy^ and Gy, is trans- 
formed into another region, denoted as Ci. In all the following, we write m to 
represent m for short and it should be noticed that 777 in this paper is not the 
Stieltjes transform as we usually denote. 

When the spiked structure (|1.1|) exists, by equation ()2.3p . this time the Stielt- 
jes transform 777 ~ rrin of F^^ satisfies 



1 1 P~M 7/„ 1 
m p 1 + 7)7 


k 
Vn V^ aiTli 

P ^ 1 + 0,777 




,^ /I p- M 7/„ J/„ ^ af 77i \ 


dm . 


\m'^ p (l + 7n)2 p ^(1 + 0^777)2 y 


'ating the same computation as before, we get: 


py-'^'-if) 


^ ^ I fi-.\ ( ^ P-M y„ Vnsr^ a 


^rii 
aimY 


p2mJc^"''™'\m'^ p {l + mf p ^ (1 + 


- ^ ^ J) ff ^ \ ^" y^- \^ {l~ai)ni \ 
p 2TTi Tc'' K m ' 1 + m p ^ (1 + 7?7)(1 + aim) J 


^ / 1 Vn , Vn Y^ 


[ 1 af 1 


j d777 


V7772 (1 + 777)2 p Z_j 


(1+777)2 (l+a,777)2 



dm 



The term 



k 



yri_sr^ (1 - ai)ni 

P p^ (1 + 777)(1 + ai77l) 
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is of order O(-), so we can take the Taylor expansion of / around the value of 

— — + 1 ^" ■ and the term 

7n l-\-ni ' 



Vn sr^ 



(1 + to)^ (1 + aimY 
is also of order O(-), this gives rise to: 



^•^'"■«-;i//'-^ 



"" -M-. 



ynm 



m 1 + m m (1 + ?7i)^ 



p Zm Jq m 1 + m p ^-^ 



i=l 
k 



(1 + ?7l)2 (1 + aiTTlY 



mdm 



p Zm Jq m 1 + m p ^-^ 



(1 - ai)ni 



I 1 2^"™ \A 



^ (1 + m)(l + fliTTl) TO (1 + ?7l) 



-0( 



1 



(2.8) 



First we consider these Ca^ii ~ I,-'- ,fc)- From the formula p.Sp . we see 
that only the last three terms contribute to the integral on the contour Ca^ (i = 
1, • • • , fc) (because for the first term, the only poles: m = and m — ~\ are not 
in the contour Cq . ) : 



n 1 

p 2Tri 



Jca, m l + m' p f-^ 

-^i fi-- + T^)tn^^^^^dm 
p Zmn Jq m 1 + m '^ — ' (1 + aiUiy 



(1 + to)2 (1 + a,m)2 



%dm 



2'Kip jQ^ 
n. 



(™+i) 



1^2 



-dm 



P 



/(0(a.)) - / (0(aO) a. 



ynflj 



/(- 



y« X yn 



/(-- 



m 1 

1 



Vn 



n 1 

p 2iTi jQ^ 

'2'^1-P Jc^^ 

-nj {(j}{a,))(a, 

p V (a,; - 1) 



(a. - 1)2 
fe 



(1 - ai)7ij 



(1 + to)(1 



ynm 



+ m p ^-^ (1 + TO,)(1 + aim) m (1 + to)' 



:)dm. 



71^(1 - a^) 



■(^ 



ynin 



TO 1 + TO (r?i + ^)ai(TO + 1) f7i (1 + to)- 

yntti \ 

l)2j- 



;)dm 



So, combining these two terms, we get the influence of the spiked part, that is, 
the integral on the contours lJi=i •••fc^ai- 



1 '' 



(2.9) 



1=1 
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So in the remaining part, wc only need to consider the integral along the 
contour Ci . Consider the second term of (|2.8p with the contour being Ci : 



p Zm Jq^ m 1 + m p ^ — ' 



k r ^ 



i=l 



(l + m)2 {l + a,m) 



mdm 



' ' n-^^ 



2'Kip 



k 



1 Mmyn M 1 Af ■^-^ riiafm 



j/„ (1 + to)2 m 2/„ m ~^ (1 + Oim) 
^G, {f)^^i /(_l + ^!^)fJL„V ,/'"'"^,, V^ (2-10) 



Combining Equations (|ZH) . (|2^ and (|2?T0| . we get: 



fc 



1 / ., 1 , Vn ., M ^r-^ riiafm 



^•'^ 2i:ipJc/ m 1 + m'^myn ^ {I + a,m)^ ' 

_^ ^ / // 1 ^ ^n ^y^ (l-aQrt, 1 y„m ^^ 

27rip Jcj ?Ti 1 + m ■^-^ (1 + ?ti)(1 + a^m) tti (1 + m)^ 

+ (1 - -)GyAf) + E -/W«^)) + 0(A) • 



«=1 P 



The proof of the theorem is complete. D 



Proof, (proof of Theorem[2]) When the a'iS satisfy |ai— 1| < y/y, from lBaik and Silverstein 



2006h . we know that the extreme eigenvalues of Sn converge to the endpoints 



of the support of the M-P distribution, i.e ay or by. So the term - X]i=i iT'if{o,i + 
^"^^ ) in Equation p.7p . which corresponds to the isolated spiked part is absent, 
and the contour C is now exactly the contour Ci which contains only the inter- 
val [j^, ^)U(^, T^] when y > 1 and [^, ^] when < y < 1. 
Taking Equation (|2.10p into consideration leads to the result. D 

3. About the contour C and Ci appearing in the main theorems and 
proofs 

As we can see in the derivations of Theorem[T]and[5J one of the most technical 
points concerns contour integrals on C and Ci . In this section, we provide more 
details on these contour integrals. 

When no spike exists, we have the equation: z = — — + -rr^- So the transform: 
z — > m{z) maps a region in the complex plan that includes the support of Gy^ 
and Gy to a region of m, which we denote as Ci. This mapping is not easy to 
visualize in the complex plan; we can consider its restriction to the real domains, 
i.e. for real z and m, as illustrated in Figure [l] with left panel corresponds to 
< y < I and right panel y > 1. Because the two situations lead to different 
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0<y<1 



y>1 




I 



-1/(1 -,/y) -1/(1 +,fy) 



-i/(i+#) -i/(i-Vy) 



Figure 1: The graph of the transform z{m) ~ 



y 

l+m ■ 



appearance of the graph, we should consider them separately. Both two panels 



have two local extrema located at (- 



and (- 



, 6) , where a and b 



are the endpoints of the support of the M-P distribution. 

The z's such that z' ( m) > are not in the support of F"^"(or Gy^), refer 
to ISilverstein and Choil (jl995l ) for further reference. First consider the case of 
< y < 1: from the left penal of Figure [H we see that the z's such that 
z'{m) > corresponds to {z > b}[J{0 < z < a}[J{z < 0}. Thus, the support 
of F^" consists of {a < z < fe}lj{~ = 0}: ^^^ changing to the variable m, 
it is equivalent to m G [y^-^, -^ 7 V ]- From the figure, we see that the point 
{m = ^71 } (equivalent to {z = 0}) is contained in this interval. As a result, 
restrict the contour Ci to the real axes, it encloses the interval 



when < y < 1. Besides, {m = —1} is the pole. For the case oi y > 1: from the 
right panel of Figure [TJ {z'{m) > 0} corresponds to {z > 6}1J{— oo < z < a}, 
thus, the support of F'^"(or Gy^J consists oi {a < z < b}, which is equivalent 
to TO € [tT^' iTV ]- Since in this situation, {z = 0} ^ {a < z < b}, so the 
value of TO = — ^ (equivalent to z = 0) should be excluded from [73-7=7 T^~7^]" 
For this reason, the contour Ci encloses [ , ~ V 1 —^) I J(^'-t: , ~V ] when w > 1. 
Besides, this time the pole is {m = 0}. 

Furthermore, suppose there exists a spiked structure as ()l.ip . so the equation 
becomes 



P 



1 



aim 



the transform: z — > m{z) maps a region of z (containing supp{F^") and 
supp{F^)) into C. This time, whether the spikes are large or small leads to 
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different situations. When ja^ — 1| > ^ (large spikes), the spike eigenval- 
ues Pi's are away fr o m th e critical values 1 ± ^/y, and it has been proven in 
Baik and Silversteinl ( 2006[ ) that the support of F^" is made of A; + 1 absolutely 



continuous components exactly: a bulk part clustering around the support of the 
M-P distribution Gy^ and k isolated small components with support [aai,bai] 
clustering around the value a^ + yai/(ai — 1) = 0(ai), i = 1, . . . , fc, respectively. 
Moreover, due to the exact separation property, F^"'^"[aai, 6aJ = "-i/Pi ex- 
actly for large p. For this reason, after mapping the support of F^" to C, it also 
consists of fc + 1 parts: Ci plus Ca^ (i = 1, • • • , fc), with each of these contour non- 
overlapping. And Cai is the contour that encloses the point of {ra = — -} (be- 
cause the k isolated parts surround around the value of z = (j){ai){i — 1, ■ • • , fc) 
and mapping these values to m is exactly rn = —^ (i = Ij • • • j^))- When 
|a,; — 1| < y^ (small spikes) . ISaik and Silversteinl ( 20061 ) shows that the extreme 



eigenvalues of S'„ cannot be separated from the support of the M-P distribution 
(they stick to the endpoints of the support: a or 6). So the contour C is now just 
Ci (without the other k components Cq. (i == 1, • • • ,k)). Moreover, since 



M 



ynsr^_ 



aim 



we see that m — ^-^ is the pole. Due to ja^ — 1| < y/y, we have always — ^ S 
\— , "'" _ , — , Z _ ] or — — e [— -I ,"'" _ , ^!-r) I \(^-T, — 1 "'" ^ 1 according to whether 
0<y<lory>l. 

We may summarize all these findings in the following Table [T] and [2] Besides, 
it should be noticed that for a given function /, there will be other poles which 
belongs to / that contribute to the integral if these poles are in the region of 



Table 1 






0<y < 1 


y>l 






^ Ml R ^ 


i 1 


region ot L-i mL[ i_^. i+yjjJ mL[ 


l+^^' i/-l''UVy_i. 


l-^^J 


Table 2 


poles in Ci < j/ < 1 


y>l 




large spikes m = —1 


m = 




small spikes m = —l,m = — — (i = !,■■■, k) 

a,; 


m = O.m = — —(i = 


--!,■■■, k) 



Last, we compare the results of Theorem[T]and Theorem[5] Seemly, Theorem[T] 
has one term more than Theorem[2] (the term - X]i=i f^ifio-i + a"-i ) ^^ missing 
in Theorem [5]) . But according to Table [21 we know that when doing contour 
integral of small spikes, the terms that correspond to the contribution of poles 
m = — — {i = !,••• jk) in Ci should be added. And it is exactly the missing 

term - X]i=i ^ifi'^i + ^""1 ) appearing in Theorem [T] if we do the calculation. 
So, Theorem [1] and Theorem [5] lead to the same result in spite of different 
expressions. Above all, we have the knowledge that whether the spikes are large 
or small has no influence on the centering parameter F^"'^"(/). 
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4. An application to the test of presence of spike eigenvalues 

In iBai et al. I ( 2009I ). a corrected likelihood ratio statistic L* is proposed to 
test the hypothesis 

Hq : T, — Ip vs. Hi : T, ^ Ip . 

They prove that under Hq, 

L*-pGyM^N{m{g),v{g)), 
where 

L* = trSn~log\Sn\-p, 

Vn 

, . log(l - y) 
Ma) = ^ , 

v{g) = -2log{l-y)-2y. 

At a significance level a (usually 0.05), the test will reject Hq when L* — 
pGy^{g) > rn{g) + ^~^{1 — a)y^v{g) where $ is the standard normal cumu- 
lative distribution function. 

However, the power function of this test remains unknown because the dis- 
tribution of L* under the general alternative hypothesis Hi is ill-defined. Let's 
consider this general test as a way to test the null hypothcse Hq above against 
an alternative hypothesis of the form: 

Hi : E has the spiked structure (|1.1|) . 

In other words, we want to test the absence against the presence of possible 
spike eigenvalues in the population covariance matrix. The general asymptotic 
expansion in Theorem [1] helps to find the power function of the test when the 
spike eigenvalues are large. However, it will lead to the same result if we use 
Theorem [2l 

More precisely, under the alternative Hi and for f{x) ^ x — logx — 1 used 
in the statistic L* , the centering term i^i'">^"(/) can be found to be 

k I- 

1 . , A/ 1 ■ 1 1 

1 + -y^n^Gi V'rijlogai - (1 )log(l-2/„) + O(^) , 

P ~l P P ~~^ Vn n^ 

thanks to the following formula 

1 . ^' M 1 

py-'"" (x) = 1 + - y n,a, + 0( — ) (4.11) 

p ^-^ p n^ 
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and 

py-"- (logx) = i Vnaoga, - 1 + (1 - — )log(l - y„) + Oi^) . (4.12) 
P ~-{ Vn n^ 

The details of derivation of these formula are given in the Appendix. 
Therefore we have obtained that under iJJ', 

It follows that the asymptotic power function of the test is 



(a) = 1-$ $"'(l-a) 



V-21og(l-2/)-22/ 



In the particular case where the spiked model has only one simple spike, i.e. 
fc = 1, ni = 1, the above power function becomes 



/3(a) = 1-$ $"'(l-a) - 



-1/-, „^ ai-1-logai 



v/-21og(l-2/)-22// ' 



which is exactly the formula (5.6) found in lOnatski et al. I ( 2011 ). Note that 



these authors have found this formula using a sophisticated tools of asymptotic 
contiguity and Le Cam's first and third lemmas. Our derivation is in a sense 
much more direct. 



Appendix A: Additional proofs of (|4.1ip and (|4.12l) 

Recall that the LRT test works only when < y < 1, so the contour C\ 
encloses the interval [j^^i TTl^] °^ '^^ ^^^ axes. And poles of {to = —1} 
and {to = — ^} (pole of the function logz) are included in this interval. In 
all the following, we only consider large spikes and do calculations according to 
Theorem m 

A.l. Proof of (|4lT|) 
We have 



IB = — 



1 / 1 i/n M M--^ riiafm 



( 1" ^i )( / T, 7TT)dm 

m 1 + 771 y„TO, -^ — ' 1 + OiTOr 



2TTip Jq_^ m 1 + 771 y,i77i ^—^ (1 + aim) 



*= -.2, 



'^T^ip Jci \m{m + 1) ^" :^ (1 + a^m)^{l + m) J 
P P ^ (l-ai)2 
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12 



^M 



"P JC^ j-[ (1 



(1 - ai)ni 



2Trip 

k 



1 



Vn-m 



+ aim){\ + m) \m (1 + m) 



\ 1 {^j / 'YYl 

» ^-^ 2 dm'' \1 + ai 

^ i=i ' 

fe 

= -E 



dm 



?n— — 1 



j=l 



^i'^iVn 



{l-a,Y 



(A.14) 



,„^^ i\/ 1 v^ / 

(EH = 1 + - > n, a, 

Combine (|A.13|) . (|A.14p and (jA.lsp . we get: 



^]+0(A). (A.15) 



F^'-^-Ca;) = 1 + - V n,a, - — + O(^) 



i=l 



A. 2. Proof of ^^ 

We have 



-1 /■ logC^j + iogl^:^;;^^),., Jk „.„|a„„,^ 



■(M-E^ 






)d'm 



Ik 9 



2mpyn Jci m 

^ A + B 



Jci m 2inpyn Jci m + 1 ~{ {^ 



dm 



+ atraY 
(A.16) 



A 



-M f ,m 
log(- 



1 



2TTipyn Jci m+1 

M f "i 

log™ • d\og{- 



) • d log m 

1 



^;.^-l 



M y, 



m.+ 1 
logm 



27ripy„ ?;„ - 1 Jc^ (m + l)(m - — ^) 



dm 



log(l-y„) , 

PVn 



(A.17) 
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1 fe 

r -"o. —)y 

2'nip Ji2^ m + 1 -^-^ (1 + airnf- 



'^Triple, "^ + 1 ^(l + a'-^2 



where 



= > (p logf ^^^^^ — )niai( -. ^)diTi 

2TTip^Jc^ ^^ m+1 ' ^l + a,7n {1 + aim)^ ' 

= C-D, (A.18) 



Imp ^-^ Jq^ m+1 1 + aiWi 

1 x~^ / ™ ^ "'_^ 

= Tr—Z^f '^jlogl p^) •dlog(l +a,m) 

27rip ^ /ci m+1 

— 1 ^ y m ^ 

= TT—y^f n,\og{l + a,m)- d\og{ ^^) 



i=l 



fe 



27rip y„-l^Jc^ (TO + l)(m-^) 



fc -, fe 



and 



= iVn,log(l-a.)-iVn,log(l + ^^), (A.19) 



D = TT—y^f log T-^ 7TT ^I^™ 

27ri|) ^ /ci m+1 (1 + aim)^ 



-y f r^^^-diog( ^) 

«P ^ Jr. 1 + ai"7 m + 1 



= yjl V/ - dm 

2Tiip{yn - 1) frt /ci (1 + a^m){m - j^){m + 1) 

Combine ((XT6)) . ((XT7| . COS]) . (|XT9|) and (|X20l) . wc get: 

AT 1 1 

(USD = -log(l-2/„) + -y]"^log(l-a^)--Vn^log(l + ^^) 

PVn p fr^ p fr^ yn - 1 

--V "', +ly^^. (A.21) 



'i=l Vn-l -^ i=l 
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Then, we consider the part p.6|) in the general formula: 

k 



2'Kip Jq^ m 1 + 771 -^-^ 1 + aim 1 + m m (1 + fT^j 

1 V^ / TO(m + l) , Oi 1 N/ 1 Unm ., 

= — y Hi (p -(:; )( -:; rT^jdm 

2'Kip ^-^ Jq^^ yn'm — m — 1 1 + aim 1 + m m (1 + ^t^J 

A; 



i--J2n.{E-F-G + H), 



2TTip{yn - 1) ^ 
where 

E^6 °'("^ + ^) , = 2m- ^""' 



Ci (l + aim)(m- ^-^) y„ + a^ - 1 ' 



Ci (to + 1)(1 + aim)(m- ^^) a.j - 1 y„ + a^ - 1 

E) ^ — = 27ri , 

r, m T 

/ U Tfl 

Collecting these four terms, we have: 

(EH = - E(^ - "^ 1 )"^ • (A-22) 

P ^ flj - 1 y„ + flj - 1 

Finally, using the know n result that Gi^^ (lo g x) = (1 — —) log(l — y„) — 1, which 
has been calculated in iBai and Silversteinl ( 20041 ) . and combine (jA.2ip . (|A.22p 



and (|2.7p . we get 

fc 

P 



py-"" (logx) = - V n, log a, - 1 + (1 - — ) log(l - y,,) + ©(^ 



i=l 
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